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We evaluate the static qqqq and qqqqq potentials in the quenched theory at /3 = 5.8 and /3 — 6.0 
on a lattice of size 16^ X 32. We compare the static potentials to the sum of two meson potentials 
for the tetraquark system and to the sum of the baryonic and mesonic potentials for the pentaquark 
state, as well as, with the confining potential obtained in the strong coupling expansion. 
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I. INTRODUCTION 

A large amount of effort has been devoted recently to 
experimental searches Q for the identification of the 
an exotic baryon state with an unusually narrow width, 
which was predicted theoretically in the chiral soliton 
model . Several studies in lattice QCD have looked for 
such a state in order to determine its mass and parity 
but no consensus has been reached yet with some groups 
reporting a bound state with mass close to the experi- 
mental value HQ and others the KN scattering state i. 
The possible existence of such a state has raised interest- 
ing questions about its structure in connection to its nar- 
row width. A number of phenomenological models have 
been put forward to explain its stability such as special 
flux tube configurations Q and diquark formation Q . 
Another exotic that might exist and has been proposed 
in the past |^ is a bound state of two quarks and two 
antiquarks. In this work we calculate the static qqqq and 
qqqqq potentials. The study of the pentaquark potential 
is particularly important for our understanding of the un- 
derlying structure of the pentaquark bound states that 
have been detected in various experiments P, 0| . 

To evaluate the static tetraquark and pentaquark po- 
tentials we construct the Wilson loops for the qqqq and 
qqqqq systems. In the strong coupling approximation 
minimization of the energy of a system of two quarks 
and two antiquarks requires that the two quarks and the 
two antiquarks are connected by the minimal length flux 
tube. The flux tubes from each of the quarks and the an- 
tiquarks can meet at one or at two points. For the geome- 
tries considered in this work the length of the flux tube 
with two Steiner points is always smaller than the config- 
uration with the flux tubes meeting at a point and there- 
fore it corresponds to the minimal flux tube length. If 
the minimal flux tube length is denoted by L^i^ then the 
confining potential in the strong coupling approximation 
is aLmin where for a we take the string tension extracted 
from the quark-antiquark potential. Besides comparing 
the static tetraquark potential to the one extracted in the 
strong coupling approximation we also compare it to the 
sum of two meson potentials. Similarly the pentaquark 
potential is compared to the potential extracted in the 
strong coupling expansion and to the sum of the bary- 
onic and mesonic potentials. Minimization of the energy 



of a pentaquark system with the condition that the flux 
tubes connecting the quarks meet at three points gives, 
as minimal length, the flux tube configuration where the 
three flux tube junctions are Steiner points. We do not 
consider here geometries resulting in the diamond flux 
tube arrangement for which the flux tubes from the 
four quarks meet at a single point where the antiquark 
is located. Obtaining results for both the tetraquark and 
pentaquark potentials enables us to look for differences in 
their behaviour which can reflect different structures in 
the tetraquark and pentaquark systems. Such informa- 
tion may lead to important insight in our understanding 
of the structure of the 8+. 



II. WILSON LOOPS 

The SU(3) Wilson loop for the tetraquark system is 
constructed by creating a gauge invariant four quark 
state at time t = 0, which is annihilated at a later time 
t as shown in Fig. ^ The two quarks are combined into 
a colour 3 and the two antiquarks into a colour 3 rep- 
resentation of SU(3). The expression for the tetraquark 
Wilson loop is given by Q 
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where the two quark lines are created at x and the two 
antiquark lines at y at time t = and annihilated at 
x' and y' at time t respectively. The staples C/(x, x',A:) 
involved in the definition of the Wilson loop are given by 



C/(x, x', k) = Pexp 



(2) 



where P is the path ordered operator and F^ denotes the 
path from x to x' for quark line k as shown in Fig. ^ 
The baryonic junction at x is joined to the anti-baryonic 
junction at y by 



f/G(x,y) = Pexp 
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at i = 0. ?7G'(y',x') is the corresponding arrow joining 
the anti-baryonic and baryonic junctions at time t. 
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The tetraquark potential is then extracted in the stan- 
dard way from the long time behaviour of the Wilson 
loop: 
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The pentaquark Wilson loop is constructed in a simi- 
lar way: The gauge invariant state evolves now two bary- 
onic and one anti-baryonic junction. The four quarks are 
grouped into two diquarks each in a colour 3 represen- 
tation of SU(3) and the two diquarks are then combined 
into a colour singlet with the remaining antiquark. The 
Wilson loop is shown in Fig. 12 and it is given by Q 
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using the same notation as that of Eq. 

III. LATTICE TECHNIQUES 

A number of improvements are needed in order to re- 
duce the noise in the measurement of the Wilson loops. 
Since the tetraquark and pentaquark potentials are larger 
than the baryonic potential they are harder to compute 
and noise reduction methods are essential. We describe 
briefly the techniques that we use in order to reduce noise 
and extract reliably the ground state: 
1. We use the multi-hit procedure replacing the tem- 
poral links by their average value 
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with SiiU) = j^TT{UA{n)FAn)) and F{n) the staple at- 
tached to the time link that is being integrated over. The 



FIG. 2: The Wilson loop for the qqqqq system. 

integration over the links in Eq. ^ is carried out analyt- 
ically It has been shown in SU(2) that replacing 
the time links by their average value in this fashion re- 
duces the error on large Wilson loops of the order of ten. 
The factor found in ref. [3 is x^* ~ 0.889^* where t is 
the time extent of the Wilson loop. For the tetraquark 
and pentaquark Wilson loops the reduction factor will 
be a;^* and x^* respectively giving an even larger noise 
reduction for the large loops. 

2. To maximize the overlap of the trial state with the 
four or five quark ground state we use APE smearing of 
the spatial links Each spatial link is replaced by 

a fat link by acting on it with the smearing operator S 
defined by 

SUj{x) =rlUj{x)+aJ2 [Uk{x)Uj{x + ak)ul{x + aj) 

(7) 

where V denotes projection onto SU(3). This is iter- 
ated n times. We consider AI different levels of smearing 
and construct an M x M correlation matrix of Wilson 
loops. We take a — 1/2, M = 3 and successive number 
of smearings ni = 0, n2 = 15 and na = 30. 

The correlation matrices C{t) for the various Wilson 
loops are analyzed using a variational method (ISj . We 
analyze the results in two different ways: 

In the first variant, which we refer to as variational 
method 1, we solve the generalized eigenvalue problem 

C{t)vk{t) ^ Xk{t)Cito)vk{t) fc = l,...,M , (8) 

taking to/a = 1 and extract the potential levels via 



aVk = lim — 
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by fitting to a constant in the range where aVk becomes 
time independent (plateau region). 

In the second variant that we call variational method 
2 we first solve the eigenvalue equation 



C{to)vk{tQ) = Xkito)vk{to) 
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and project the Wilson correlation matrices to the space 
corresponding to the N largest eigenvalues 

C^^it)^iv,ito),C{t)v,ito)) (11) 

with N < M. We then solve the generalized eigenvalue 
equation 

C''{t)vk{t) = Xk{t)C''{to)vk{t) (12) 

in the truncated space. 

In addition the potential is extracted by considering 
the Wilson loop with the largest number of smearings 
and fit to the ratio 

V= lim- \n[W{t)/W{t - 1)] (13) 

t — 'OO 

in the plateau region. 

The noise reduction techniques described above were 
shown to yield accurate enough results for the baryonic 
potential for inter-quark distances of 0(1.5) fm [T^. 

IV. RESULTS 

All the computations were carried out on a lattice 
of size 16'^ X 32 at /3 values 5.8 and 6.0 using 200 and 
220 configurations respectively available at the NERSC 
archive |17l |. 

We compute the static qq potential on the same con- 
figurations used for the evaluation of the tetraquark and 
pentaquark potentials. Fitting the qq potential to the 
Cornell Ansatz 

V,g{r) =Vo-^ + ar (14) 

we extract the parameters Vq, a and string tension a. 
We give the values obtained in Table Using the 
value of a^CT given in Table HI and the physical value of 
y/a — 440 MeV known from Regge theory enables us 
to fix the lattice spacing a. We obtain a = 0.10 fm at 
f3 = 6.0 and a = 0.15 fm at /3 = 5.8. Since in this 
work we compare the tetraquark potential to the sum 
of two meson potentials and the pentaquark potential 
to the sum of the corresponding mesonic and baryonic 
potentials we need to compute all the potentials with 
the quarks at the same locations. For the mesonic po- 
tential needed in the comparison of the tetraquark sys- 
tem we use lattice data avoiding any parametrizations. 
For the baryonic potential needed in the comparison of 
the pentaquark potential we do not have lattice data 
and therefore we need a parametrization. Two possi- 
ble Ansatze have been discussed in the literature for the 
baryonic potential, the Y- Ansatz 18] and the sum of two- 
body potentials or A- Ansatz. The latter was derived 
using center vortices (19j but with an erroneous assump- 
tion, which after correction was shown to support the 
Y- Ansatz 20|. Nevertheless within lattice QCD one can 
check the two Ansatze and determine whether the con- 
fining potential is closer to a sum of two-body potentials 



TABLE I: The parameters of the static qq potential in lattice 
units at /3 = 5.8 and /3 = 6.0. 
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0.637(5) 


0.255(5) 


0.050(1) 


5.8 


0.636(11) 


0.248(11) 


0.105(2) 



or to the Y- Ansatz as a function of the distance between 
the quarks. The results of the two Ansatze differ at the 
most by 15% and have been compared to lattice results 
in refs. [S IH IH IH IM El ■ The general consensus 
from these studies is that at large distances the bary- 
onic potential approaches the Y- Ansatz. Therefore here 
it suffices to use the Y- Ansatz to parametrize the bary- 
onic potential. As in rcf. 16] we use the values for Vq, a 
and a extracted from fitting the qq potential via Eq. H14|) 
and given in Table This means that there are no ad- 
justable parameters in the comparison of the pentaquark 
potential to the sum of the mesonic and baryonic po- 
tential. In addition both the tetraquark and pentaquark 
potentials are compared to the potential derived in the 
strong coupling expansion. We use the parametrization 
given by 




(15) 

where again for Vq , a and a we use the values given in 
Table H The factor Uq in front of the Coulomb term is 
one when a quark interacts with an antiquark and one 
half when the interaction is between quarks or between 
antiquarks as obtained from the one-gluon exchange ap- 
proximation. 
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FIG. 3: The tetraquark minimal length. The angles at the 
junctions are 120'\ 

For a general tetraquark configuration the minimal 
length having two Steiner points is shown in Fig. Inland, 
in general, is non-planar. For the planar geometry con- 
sidered in this work the minimal length can be easily 
computed and it is given by 

^min = + V3R2,R2 + VSRi^ , (16) 

where the distances i?i and i?2 are defined in Fig. 01 For 
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FIG. 4: The geometry used for the computation of the 
tetraquark static potential where we have taken the bary- 
onic junction at the origin and the anti-baryonic junction 
at y = (_Ri,0, 0). The quarks are located at positions 
ri = (0, -R2/2, 0), r2 = (0, — ii2/2,0) and the antiquarks at 
rg = (J?i,i?2/2,0) and = (iii, -i?2/2, 0). 




whereri2 = |ri-r2|,r23 = |r2-rs3 1, and rai = |rs3-ri|. 
The minimal length joining the three points ri, r2 and 
is equal to B. If any interior angle of the triangle 
rir2i"S3 is greater than 120" then the minimal length is 
given by ri2 + r23 + rai — max(ri2 , r23 , r3i ) . This is taken 
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FIG. 5: The pentaquark minimal length. The angles at the 
junctions are all 120°. 

a general pentaquark configuration the minimal length 
with three junctions is non planar and has the geometry 
shown in Fig. [S] with the three junctions being Steiner 
points. We compute these three Steiner points numer- 
ically as follows: We start by an initial guess for the 
Steiner point S3. Having an initial guess for S'3 enables 
us to evaluate the two other Steiner points analytically 
in terms of S'3 following ref . J8j : 

- ife 

- rsi I = («, j, k cyclic) 

^ — '^12 + ^23 + ^31 



FIG. 6: Plateau values for the tetraquark potential. Top for 
the case where Ri/a — 2 and R2/a — 8; middle for Ri/a — 4 
and R2/a = 4; bottom for Ri/a — 8 and i?2/a = 2. The 
data obtained using variational method 1 are denoted by the 
circles, using variational method 2 by the filled triangles and 
using Eq. by the crosses. The circles and filled triangles 
are shifted horizontally for clarity. 



into account in our analytic determination of the Steiner 
points. An analogous expression holds for 6*2. Having 
a value for Si and S2 and knowing rs we can evaluate 
analytically the new value for S'3. This procedure is it- 
erated until convergence is reached after a small number 
(of order 10) of iterations. We use this method to com- 
pute for the geometry studied in this work. This 
approach can also be applied to find the minimal length 
for a general tetraquark geometry. 



We discuss first the results obtained for the tetraquark 
system. As already mentioned the geometry that we use 



for the computation of the tetraquark potential is planar 
and is shown in Fig. ^ When R2 is small this configu- 
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TABLE II: Determination of plateau values for the tetraquark potential at /3 = 6.0. The first column gives the diquark distance, 
ill, and for each of the internal diquark distances, R2, we give, the the initial and final time ti/a and tf/a used in fitting the 
plateau, the value of the potential and = x'^/d.o.f . We denote with an asterisk the value that we choose. All quantities are 
given in lattice units. 
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FIG. 7: The tetraquark static potential at /? = 6.0 (crosses) 
and l3 — 5.8 (filled triangles) for R2 = 0.6 fm as a function 
of the minimal length in physical units. We have applied a 
constant shift to the data at /3 = 5.8. The jackknife errors 
are smaller than the size of the symbols. 



ration allows diquark formation. We take the distance 
between the two quarks and the two antiquarks, i?2, to 
be equal. We will refer to R2 as the internal diquark 
distance. The potential is computed as a function of the 



distance between the two diquarks, and the internal 
diquark distance i?2- In Fig. El we show the plateaus for 
representative values of the distances i?i and i?2- On 
these figures we display the values extracted using the 
two variational methods as well as the one extracted from 
the ratio of Wilson loops with the largest number of APE 
smearings as given in Eq. H13|) . As can be scon from these 
figures the two variational methods yield consistent re- 
sults of similar quality. The ratio of the Wilson loops 
on the other hand yields better plateaus. This is true 
for all the distances that we have performed this analy- 
sis and we will therefore use this method to extract the 
potential. We have always checked that the variational 
analysis gave results consistent with those extracted us- 
ing Eq. H13|l. The errors shown on all our figures are 
jackknife errors. 

We check scaling by comparing results at /3 — 5.8 and 
/3 — 6.0. In Fig.[7|we show the static tetraquark poten- 
tial for the same physical internal diquark distance R2 as 
a function of the minimal length. We apply a constant 
shift to the data at /3 = 5.8 given by twice the differ- 
ence in the constants Vq at the two /? values. As can 
be seen the results, to a good approximation, fall on a 
universal line showing reasonable scaling behaviour. As 
in all our figures the errors shown on this figure are the 
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w(t-i) ^'^'^ tetraquark static 
potential versus time in lattice units. Upper graph for 
Ri = 0.2 fm and R2 = 0.6 im ed, (3 — 6.0. Lower graph 
for Ri — 1.2 fm and R2/a — 0.6 fm at /3 = 5.8. Horizontal 
lines show the value obtained from fitting the lattice data to a 
constant with the solid line showing the value that we choose. 
The initial time used in the fit is shown by the position where 
the line begins. The dotted curve shows the result of fitting 
to Eq. [H 
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FIG. 9: The tetraquark static potential ai f3 — 6.0 versus 
the minimal length. Data obtained for R2 =0.2, 0.4, 0.6 and 
0.8 fm are denoted with crosses, filled triangles, filled squares 
and filled circles respectively. 
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statistical errors obtained by jackknife analysis. To un- 
derstand why the statistical errors are small even at the 
larger distances we plot in Fig. |S1 for two representative 
cases, the ratio Vcs{t) — — log , which for large t 

gives the ground state potential V^q we are interested in. 
The two cases that we consider are those that have the 
smallest and largest values of Lmin in Fig. [7| The first 
corresponds to Ri — 0.2 fm a.t f3 — 6.0 and the second 
to Ri = 1.2 fm at /3 = 5.8 with R2 = 0.6 fm in both 
cases. To fit to a constant we must search for a plateau 
of Vee{t) and make sure that changing the initial fitting 
range does not produce a value outside the statistical er- 
rors. In our analysis we take the plateau value that gives, 
for the earliest initial time range, = x^/d.o.f. ~ 1 
since this criterion ensures a good fit and gives a result 
with the smallest statistical error. When i?i = 0.2 fm the 
results are very accurate and clearly show contributions 
from excited states with Vce{t) showing a plateau only 
for t/a > 7 . Fitting within the plateau range we extract 
the values given in Table UTI where we check consistency 
by changing the initial fit range. Using our criterion we 
choose the plateau value with = 0.8 denoted by the as- 
terisk in Table Ull At these small distances the accuracy 
of the data allows, in addition, a fit that includes the first 



FIG. 10: The tetraquark static potential versus Ri for 
R2 =0.2 fm (crosses), 0.4 fm (filled triangles), 0.6 fm (filled 
squares) and 0.8 fm (filled circles). The solid line is twice the 
qq potential. 



excited state. This enables us to extend the time range of 
the fit and check that including the first excited state re- 
produces the plateau value obtained for the ground state 
from fitting to a constant at larger times. In this case 
the time-dependence of Vcs{t) is given by 



V,s{t) = Viq - log 



1 



-dVt 



1 + ce-''^(*-i) 



(18) 



where dV is the gap between the ground state and the 
first excited state. Fitting from ti/a = 3 we obtain 
aViq — 1.195(1), which agrees with the value extracted 
by fitting to a constant within the range of the plateau. 
When we consider the second case at the largest value of 
Lniin in Fig.[71we see from Fig.|Slthat 14ff (0 has very large 
errors when t/a > 5. This means that the first accurate 
points used in the fit will determine the value and error 
of the extracted constant and that checking for excited 
state contributions is no longer possible. In this partic- 
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FIG. 11: The tetraquark static potential at /9 = 5.8 
for i?2 = 0.3 fm compared with V^'^ (sohd line) and 
with 2Vqg(i?i) (dashed line). 
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(fm) 

FIG. 12: The tetraquark static potential at /9 = 6.0 
for i?2 = 0.8 fm compared with V^^^ (solid line) and 
with 2Vqq{R\) (dashed line). 



ular example fitting to a constant from tija = 3 gives 
< 1.0 with an error on the extracted value which is 
small since it is basically determined from the two first 
accurate points. Changing the initial fit range to t^/a = 5 
lowers the plateau value as shown in the figure but car- 
ries an error similar to that of the data point at t/a = 5 
making it consistent with the value extracted from the 
fit that used ti/a = 3. Within our fit criterion we choose 
the value obtained from the fit with ti/a = 3, which car- 
ries a small error explaining the small error of the data 
point at the largest value Lmin in Fig.[7| We note that for 
these large interquark distances changing the final time 
range will not affect the value and error of the extracted 
potentials since the Wilson loops become very noisy at 
the larger times. In this particular example whether we 
take the upper fit range tobei//a = 7ori//a = 12 does 
not affect the value and error resulting from the fit. As 
can be seen from Table Ull in general we choose the larger 
available value of tf /a but for the large distances we re- 
duce the upper fit range when the data points become 



q 




q 



FIG. 13: The geometry used for the computation of the 
pentaquark static potential where we have taken the baryonic 
junctions at x = (iii, 0, 0) and at z — (0, Ri, 0) and the anti- 
baryonic junction at y = (0,0,0). The quarks are located 
at positions ri = (i?i, 0, i?2/2), r2 = (7?i, 0, — 7i2/2), rs = 
(0, i?2/2) and r4 — {0, Ri, —R2/2) and the antiquark at 
rs = (0,0, i?i). 

too noisy. Since our accurate data points are limited 
to smaller times the plateau values that we extract may 
have excited state contributions, which can not be esti- 
mated since small shifts in the data points are not visible 
due to their large errors. However the fact that the data 
at the larger values of Lmin fall on the same curve as the 
better determined points is an indication that this error 
is reasonably small. In Table HTl we give the dependence 
of the plateau values on the fitting range and our choice 
of the value plotted in the figures. The table is done for 
(3 — 6.0 since this is the case that we discuss in more de- 
tail in the manuscript. The analysis at /? = 5.8 is carried 
out in the same way. 

In Fig. we show the results for the tetraquark poten- 
tial at /? = 6.0 as a function of the minimal length. As 
can be seen we have four different sets of curves for each 
internal diquark separation i?2=0.2, 0.4, 0.6 and 0.8 fm. 
When the distance between the diquarks becomes greater 
than the internal diquark distance the four sets of curves 
converge to the same line. This means that for i?i > R2 
the tetraquark potential becomes approximately only a 
function of On the other hand when Ri < R2 and 

for i?i ^ 0.4 fm the potential is independent of i?2- This 
behaviour is better seen in Fig. 1101 where we show the 
tetraquark potential as a function of i?i for the four dif- 
ferent values of the internal diquark distance R2: When 
Ri = 0.2 fm the potentials coincide for all four different 
values of i?2- When i?i = 0.3 fm and, to a good approx- 
imation, also when i?i = 0.4 fm the potentials coincide 
only for R2 = 0.4, 0.6, 0.8 fm. In fact, for i?i = 0.2 fm 
and 0.3 fm and R2 > Ri, the tetraquark potential is 
given by the sum of the two meson potentials as can be 
seen by the agreement of these results with the line on 
the figure showing 2Vqq{Ri). This means that for these 
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FIG. 14: Plateau values for the pentaquark potential. Top for 
the case where iii/a = 1 and Kija — 8; middle for Ri/a = 4 
and R2/a = 2; bottom for Ri/a — 8 and R'i/a = 2. The rest 
of the notation is as in Fig. |S] 
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FIG. 15: The pentaquark static potential at /3 = 6.0 (crosses) 
and /3 = 5.8 (filled triangles) for 7?2 = 0.6 fm as a function 
of the minimal length. The data at /3 = 5.8 are shifted by a 
constant. 



geometries we have a system of two mesons rather than a 
genuine four quark bound state. On the other hand when 
the distance between the diquarks is larger than their in- 
ternal diquark distance the potential is lower than the 
sum of the two qq potentials. In Fig. ^] we show lattice 
data for (3 = 5.8 for the smallest internal diquark sepa- 
ration, namely for R2 = 0.3 fm, and compare them with 
and 2Vqq{Ri). We have chosen j3 = 5.8 in order 
to reach larger diquark separations. The tetraquark po- 
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FIG. 16: The pentaquark static potential at /3 = 6.0 versus 
the minimal length. The notation is the same as that used in 
Fig. El 



tential starts as a sum of the two meson potentials and 
then crosses over to approach V^l^. Although V^l^ has 
a larger slope as compared to the slope of the tetraquark 
confining potential it approximates best the lattice data 
for values of Ri in the range of about one 0.5 to 1 fm. 
In Fig. 1121 we perform the same comparison but for in- 
ternal diquark distance equal to 0.8 fm. In this case 
the tetraquark potential becomes larger than the sum of 
the two meson potentials approaching V^f^ from below. 
For the two largest distances the data become too noisy 
and better statistics are required to determine more ac- 
curately the long distance behaviour. 

We perform a similar analysis for the pentaquark sys- 
tem. The geometry used for the pentaquark potential is 
shown in Fig. ^| We place the antiquark on the z-axis 
distance i?i from the origin. The two pairs of quarks are 
placed at distances (i?i, 0, ±i?2/2) and (0, ±i?2/2) so 
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FIG. 17: The pentaquark static potential versus Ri for 
i?2 =0.2 fm (crosses), 0.4 fm (filled triangles), 0.6 fm (filled 
squares) and 0.8 fm (filled circles). The baryonic potential 
shifted by a constant is shown by the open squares. 
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TABLE III: Determination of plateau values for the pen- 
taquark potential ai f3 — 6.0. The notation is the same as 
that of Table ini 



Pentaquark potential 



Ri/a 


i?2/a = 2 


i?2/a = 4 


Ri/a = 6 


Ri/a = 8 




ti/a 


tf/a 


aVsq 




ti/a 


tf/a 


aVsq 




ti/a 


tf/a 


aVsq 


x' 


ti/a 


tf/a 


aVsq 


x' 




11 


12 


1.369(1) 


2.8 


7 


11 


1.616(1) 


8.5 


5 
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1.885(2) 


2.9 
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2.111(2) 
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"This value was extracted using Eq. 1181 



that when R2 = this configuration reduces to the ge- 
ometry that we used for the baryonic potential [T^ . The 
quality of the plateaus is seen in Fig. ^1 for representa- 



tive values of Ri and R2 . Again the ratio of Wilson loops 
with the largest number of APE smearings produces the 
best plateaus. 



r 



In Fig. 1151 we compare the results for the pentaquark 
potential at (3 = 5.8 and 6.0 for the same physical in- 
ternal diquark separation R2. The data at (3 = 5.8 are 
shifted by 5/2 the difference in the values of Vq at the two 
f3 values. As in the case of the tetraquark potential, the 
results for the pentaquark potential approximately fall on 
the same curve indicating good scaling. The comments 
made for the errors on the results obtained aX j3 — 5.8 
at the larger distances when discussing the scaling of the 
tetraquark potential apply also here. Table UTTl gives de- 
tails of our fitting procedure at /3 = 6.0 and provides an 
indication of the systematic errors involved in changing 
the fitting range. 

In Fig. 1161 we plot the pentaquark potential evaluated 



at /? = 6.0 versus L^^:^^. Although not as clearly seen as 
in the case of the tetraquark system there are four sets of 
curves corresponding to the four different values of i?2- 
For large values of i?i they converge to the same curve 
which again suggests that the potential asymptotically 
depends only on the minimal length at least to a first 
approximation. The four sets of curves are more clearly 
seen in Fig. 1 171 where the pentaquark potential is plotted 
as a function of the diquark separation i?i for the four 
values of the internal diquark distance R2 = 0.2,0.4,0.6 
and 0.8 fm. On the same figure we also display our lattice 
results for the baryonic potential shifted by a constant 
showing that the linear dependence of the pentaquark 
potential for small separations R2 is approximately the 
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same as that of the baryonic potentiaL It should be noted 
that, for the geometry used in the evaluation of the pen- 
taquark potential, the sum of the baryonic and mesonic 
potentials depends on both Ri and R2 and therefore we 
can no longer draw a universal curve for all separations 
i?2 as we did for the case of the tetraquark potential. In- 
stead we compare the pentaquark potential to V^?^ and 
the sum of the baryonic and mesonic potentials in Figs. 1181 
and El for two extreme cases: In Fig. E| we show re- 
sults for the smallest possible internal diquark distance, 
namely for R2 = 0.3 fm, and consider (3 = 5.8 so that we 
can reach larger physical distances in Ri . In this case the 
potential approaches V^^^ for Ri > 0.4 fm. This means 
that when Ri becomes larger than the internal diquark 
distance the qqqqq system is well described by the mini- 
mal flux connecting the quarks. This genuine pentaquark 
state has static energy which is lower than the sum of the 
baryonic plus the mesonic potential. In Fig. 1191 we show 
results for the largest possible internal diquark separa- 
tion, namely for R2 = 0.8 fm, and consider P = 6.0 since 
the results for (3 — 5.8 become too noisy to extract any 
useful information. For this geometry the quarks in the 
diquarks are always a distance larger or equal to Ri. In 
this case the results are well described by the sum of the 
baryonic and mesonic potentials and only for distances 
larger than 0.6 fm they begin to approach V^^^ from be- 
low. More statistics are required to reduce the errors in 
order to draw a definite conclusion for the large distance 
behaviour in this case. 



CONCLUSIONS 
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FIG. 18: The pentaquark static potential at /? = 5.8 
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We have calculated the static tetraquark and pen- 
taquark potentials selecting geometries that are moti- 
vated by the diquark picture. Using multi-hit for the 
time links and APE smearing for the spatial links we 
obtain results for inter-quark distances of the order of 
1 fm. Comparing the results obtained at /? = 5.8 and 
6.0 we show that the potentials have reasonable scaling 
properties. 

The main conclusions regarding the tetraquark po- 
tential are: 1) When the two diquarks are closer than 
~ 0.5 fm and with internal diquark distances larger than 
their separation (i?2 > Ri) then the qqqq system breaks, 
as expected, into two mesons and its static potential is 
approximately given by the sum of the two meson po- 
tentials. The geometry used for the tetraquark poten- 
tial makes it easy to see when the ground state is a two 
meson state since then the static potential becomes in- 
dependent of the internal diquark distance R^. 2) For 
distances between diquarks larger than the internal di- 
quark distance the static potential is approximately de- 
scribed by l^^fn' which has a confining part proportional 
to the minimal length fiux tube joining the quarks. In the 
parametrization of V^^^ we used values for Vo,a and a 
extracted from the static qq potential and therefore there 
are no adjustable parameters. The linear dependence of 



FIG. 19: The pentaquark static potential at /3 = 6.0 
for 7?2 = 0.8 fm compared with (solid line) and 
with the sum of the baryonic and mesonic poten- 
tials (dashed line). 



the tetraquark potential would require a smaller value of 
a as compared to that extracted from the qq potential. 
For internal diquark separations between about 0.5 and 
1 fm the potential is larger than the sum of the two meson 
potentials and approaches T^J^'j from below. 

The pentaquark potential shows similar behaviour: For 
diquark separations R\ small as compared to R2 the po- 
tential is well approximated by the sum of the mesonic 
and baryonic potentials. For values of R\ larger than 
i?2 it approaches However in this case I^^f^ P'"°" 

vides a better description to the data than it did for 
the tetraquark case. For our geometry we were able 
to confirm this behaviour for internal diquark distances 

i?2 ^ 0.5 fm. For larger values of i?2 the data indi- 
cate that the pentaquark potential also approaches I^^j^ 
albeit with large statistical errors that need to be re- 
duced before one can make a definite statement. When 
i?2 ~ 0.5 fm the pentaquark potential approaches V^f,j 

from above. When i?2 ^ 0.5 fm the pentaquark poten- 
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tial approaches VJ^?^ from below at least for distances Ri 
up to ^ 1 fm studied in this work. Although the exact 
values for the distances i?i and i?2 where this behaviour 
is observed may depend on the geometry, the general be- 
haviour should not be affected. 

In summary, these observations suggest that, when 
the separation between diquarks is larger than the in- 
ternal diquark distances with the latter obtaining values 
up to about half a fermi, the tetraquark and pentaquark 



systems behave as multiquark bound states rather than 
break into mesons or baryons. Although this behaviour 
is determined using a particular geometry it should hold 
for general geometries that allow diquark formation. 
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Note Added: After the submission of this work we 
notice a preprint on the static potential on hep-lat ,26(| . 
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